Critical and non-critical jamming of frictional grains 
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We probe the nature of the jamming transition of frictional granular media by studying their 
vibrational properties as a function of the applied pressure p and friction coefficient ^. The density 
of vibrational states exhibits a crossover from a plateau at frequencies lj > ijj*{p,fj,) to a linear 
growth for u < u}*(jp,fi). We show that ui* is proportional to Az, the excess number of contacts 
per grains relative to the minimally allowed, isostatic value. For zero and infinitely large friction, 
typical packings at the jamming threshold have /S.z — > 0, and then exhibit critical scaling. We study 
the nature of the soft modes in these two limits, and find that the ratio of elastic moduli is governed 
by the distance from isostaticity. 

PACS numbers: 45.70.-n, 46.65.+g, 83.80.Fg 



Granular media, such as sand, are conglomerates of 
dissipative, athermal particles that interact through re- 
pulsive and frictional contact forces. When no exter- 
nal energy is supplied, these materials jam into a dis- 
ordered configuration under the action of even a small 
confining pressure l]. In recent years, much new insight 
has been amassed for the jamming transition of models 
of deformable, spherical, athermal, frictionless particles 
in the absence of gravity and shear 2]. The beauty of 
such systems is that they allow for a precise study of 
the jamming transition that occurs when the pressure p 
approaches zero (or, geometrically, when the particle de- 
formations vanish). At this jamming point "J" and for 
large systems, the contact number 3] equals the so-called 
isostatic value ^ (see below) , while the packing density 
(f>j equals random close packing 0, 0|. Moreover, for 
compressed systems away from the jamming point, the 
pressure p, the excess contact number Az = z{p) — zP^ 
and the excess density — (f) — (f)j are related by power- 
law scaling relations — any of the parameters p, Az and 
A(f) is sufficient to characterize the distance to jamming. 

Isostatic solids are marginal solids — as soon as con- 
tacts are broken, extended "floppy modes" come into play 
13 . Approaching this marginal limit in frictionless pack- 
ings as p — > 0, the density of vibrational states (DOS) 
at low frequencies is strongly enhanced — the DOS has 
been shown to become essentially constant up to some 
low- frequency crossover scale uj*, below which the con- 
tinuum scaling - w'^-i is recovered H H H H S H [HI . 
For small pressures, ui* vanishes ~ Az. This signals the 
occurrence of a critical lengthscale, when translated into 
a length via the speed of sound, below which the material 
deviates from a bulk solid (9|i • The jamming transition for 
frictionless packings thus resembles a critical transition. 

In this paper we address the question whether an anal- 
ogous critical scenario occurs near the jamming transition 
at p = of frictional packings. The Coulomb friction law 
states that when two grains are pressed together with a 
normal force F"^, the contact can support any tangential 



friction force with F^ < /iF", where /x is the friction 
coefficient. In typical packings, essentially none of these 
tangential forces is at the Coulomb threshold F* = //F" 
Hi ^1 • A crucial feature of these packings of frictional 
particles is that, for p 0, they span a range of pack- 
ing densities and have a nonunique contact number zj (/i) , 
which typically is larger than the frictional isostatic value 
z.'^^^ = d-H 1 (see below). So two questions 

arise: do frictional systems ever experience a "critical" 
jamming transition in the sense that lu* vanishes when 
p — > 0? What is the nature of the relations between p, 
fi, UJ* , excess contact number z{p,p) — zj{fj,) and excess 
density? 

The results which we present below give convincing ev- 
idence that jamming of frictional grains should be seen as 
a two-step process. The first step is the selection of z for 
fixed /i and a given numerical procedure (see Fig. la); this 
has been studied before [ij, [ij, |l3 . Our focus here 
is on the second step, the fact that the critical frequency 
UJ* of the DOS of vibrations of infinitesimal amplitude 
is proportional to the distance to the frictional isostatic 
point Az := z{fi,p) — Zj'^^. The crucial point is that zj(/i) 
and Zj'^Q in general differ. In particular for small values 
of n, the contact number saturates at a value substan- 
tially above the isostatic limit, uj* saturates at a finite 
value and the system remains far from criticality. For in- 
creasing values of /i, however, zj{p) approaches z^^^, and 
thus for large friction values, uj* evidences an increas- 
ingly large scale near the jamming point. The vanishing 
of UJ* has its origin in the emergence of floppy modes at 
the isostatic point. We show that similarly to frictionless 
systems, both uj* and the ratio of shear to compression 
modulus G/K scale as Az. In short, the distance to 
isostaticity, which is well-defined, governs the scaling of 
both frictional and frictionless systems, providing a uni- 
fied picture of jamming of weakly compressible particles. 

Let us, before presenting our results, recapitulate the 
well-known counting arguments for the contact number 
in the limit p — > for dimension d Since the de- 
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FIG. 1: (a) Average contact number z as a function of p^^^ for various ji as indicated, (b-e) Vibrational DOS for granular 
packings for friction coefficients as indicated, and for pressures approximately 5 x 10~^, 5 x 10~^, 5 x 10"'*, 4 x 10"'^ and 3 x 10~^. 
For decreasing p the DOS becomes steeper for small uj, and the crossover frequency uj* , indicated in (e), decreases with p. The 
packing with = 0^ is obtained by first making a frictionless packing and then turning on the tangential frictional forces in 
the DOS calculation. As noted in the text, all frequencies are scaled by a factor p^^^ ■ 



formation of the spheres vanishes in the limit p — Q, all 
particles in contact are at a prescribed distance, which 
gives zN/2 constraints on the dN particle coordinates, 
leading to z < 2d. For the frictionless case, the zN/2 
normal contact forces are constrained by Nd force bal- 
ance equations — hence only for z > 2d can we generi- 
cally find a set of balancing forces [3. Taken together, 
this yields z^2d—: as p 0: at the jamming tran- 
sition, packings of frictionless spheres are isostatic. For 
frictional packings, there are zdN/2 contact force com- 
ponents constrained by dN force and d{d — \)N /2 torque 
balance equations — thus z > rf + 1, with zt^^^ = d + 1 
the isostatic value. Hence, at the jamming transition, 
frictional spheres do not have to become isostatic but 
can attain contact numbers between z^^^ = d -I- 1 and 2d. 
While it is not well understood what selects the contact 
number of a frictional packing at J, simulations for discs 
in 2 dimensions show that in practice zj{p,) is a decreas- 
in g fu nction of /Lt, ranging from 4 at small /i to 3 for large 
51 El El; see also Fig.HJi. 

Procedure — Our numerical systems are 2D packings 
of 1000 polydisperse spheres that interact through 3d 
Hcrtz-Mindlin forces |l7| , contained in square boxes with 
periodic boundary conditions. We set the Young mod- 
ulus of the spheres E* = 1, which becomes the pres- 
sure unit, and set the Poisson ratio to zero. Our unit of 
length is the average grain diameter, the unit of mass is 
set by asserting that the grain material has unit density 
and the unit of time follows from the speed of sound of 
pressure waves inside the grains E3- The packings are 
constructed by cooling while slowly inflating the particle 
radii in the presence of a linear damping force, until the 
required pressure is obtained. For each value of and p, 
20 realizations are constructed (occasional runs with 100 
realizations did not improve accuracy). 

Once a packing is made, the additional damping force 
is switched off and the dynamical matrix is obtained by 
linearizing the equations for small amplitude motions, 
which include both rotations and translations. It is im- 
portant to realize the special role of the friction: if the 
density of contacts that precisely satisfy = is 



negligible, the Coulomb condition F' < /j,F" only plays a 
crucial role during the preparation of a packing. We will 
assume that this is the case, and come back to this subtle 
point later. Under these assumptions, and for arbitrar- 
ily small amplitude vibrations, the Coulomb condition is 
automatically obeyed and the value of /i does not play a 
role anymore in analyzing the vibrational modes. More- 
over, the changes in F* are then non-dissipative and the 
eigenmodes of the dynamical matrix are undamped. In 
this picture, the main role of the value of the friction 
coefficient is in tuning the contact number. 

We analyze the density of vibrational states (DOS) of 
the packings thus obtained. Since for Hertzian forces 
the effective spring constants scale with the overlap 6 as 
dF"/d(5~(5^/^~p^/'^ jl3], all frequencies will have a triv- 
ial p^^^ dependence. To facilitate comparison with data 
on frictionless spheres with one-sided harmonic springs 
II 13 ! we report our results in terms of scaled frequencies 
in which this p^^^ dependence has been taken out. 

Variation of z — Anticipating the crucial role of 
the contact number, we start by presenting z(/i,p) for 
our packings. Fig. la confirms the earlier observations 
iil2;i.l3i 14] that the effective value of zj(/i) =: z(/i,p ^ 0) 
varies from about 4 to about 3 when /i is increased. More- 
over, the excess number of contacts z{p,,p) — zj(/i) varies 
with pressure as p^^^ for all values of /i. 

DOS — Figs, n^-e show our results for the DOS for 
various values of fi. For the frictionless case shown 
in Fig. db), we recover the gradual development of a 
plateau in the density of states as the pressure is de- 
creased HQ- For this case z z? and the crossover 
frequency u>* scales as Az = z — zP^ tllalsl (see below). 
However, as Figs.^c-d) illustrate, as soon as the tangen- 
tial frictional forces are turned on, this enhancement of 
the DOS at low frequencies largely disappears, because 
the frictionless "floppy modes" are killed. This point is 
demonstrated most dramatically in Fig.^c), where the 
underlying packing has been generated for zero friction, 
and the friction is only "switched on" when calculating 
the DOS — this represents the limit of vanishing small 
but nonzero friction, for which the DOS is seen to be 
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FIG. 2: (a) uj* for /i = 10 as function of uiovcviap — the regime 
(1 < tJoveriap < 3) corresponds to the crossover regime in the 
DOS that we focus on here (see text and panel (b)). (b) The 
rescaled DOS for ^ = 10 exhibit good data collapse in the 
crossover regime. Here 20 rescaled DOS are shown with p 
rangmg from 9xl0~^ to 3x10"^ 



very far from critical. By increasing the friction coeffi- 
cient, the development of a plateau and the scaling of the 
crossover progressively reappear (Fig. H^e)). The intu- 
itive picture that emerges is that with increasing friction, 
granulates at the jamming point approach criticality. 

In order to back this up quantitatively, we perform a 
scaling analysis of the low-frequency behavior of these 
DOS. To avoid binning problems, we work with the in- 
tegrated density of states I{uj) — duj'DOS{Lj'). The 
critical frequencies are then obtained by requiring that 
the rescaled integrated DOS, {uj*)^^I{oj/u!*), collapse. 
Such collapse is never perfect, in particular since not 
all DOS have precisely the same "shape" (Fig. 1). We 
vary the value of uJovcviap ■= uj/u* where we require the 
rescaled integrated DOS to overlap — as Fig. illus- 
trates, this yields precise values for lo* as function of 
(^overlap- Restricting ourselves to the crossover regime 
(1 < cjovoriap < 3), wc obtain by this procedure both an 
estimate of uj* and of its errorbar. As Fig. illustrates, 
when rescaled with these estimated values of uj*, the col- 
lapse of the DOS in the crossover regime is convincing. 

Scaling of uj* — The first main result of this paper 
is shown in Fig. |21 uj* does not scale in a simple way 
with p, but the data for all fi and p collapse onto a 
single curve when plotted as function of Az = z — Zj^^ 
(Fig. |3t,c). Moreover, uj* ~ Az — the plot of z — 3 
versus p shown in Fig. is essentially equivalent to the 
plot of UJ* vs p. In other words, packings with Az 1 
have many low-lying vibration modes and correspond- 
ingly a large enhancement (plateau) in the DOS. The 
dominant quantity governing the behavior of frictional 
granular media is the distance from the frictional "crit- 
ical point" z = Zj'^Q — d + 1 — 3. This distance can be 
characterized most conveniently by Az — put in these 
terms, the scaling of lj* for frictional media is very similar 
to the scaling for frictionless media shown in Fig.|Sld,e). 

Scaling of Elastic Moduli — The contact number for 
isostatic systems reaches the minimum needed to remain 
stable — hence additional broken bonds then generate 
global zero-energy displacement modes, so-called floppy 
modes For frictionless systems, the excess of soft 

modes and development of a plateau in the DOS for 
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FIG. 3: (a) uj* as function of pressure p for a range of friction 
coefficients — errorbars are similar or smaller than sym- 
bol sizes (b) Deviation from isostaticity for the same range 
of parameters, (c) uj* scales linearly with the distance to 
isostaticity for frictional packings, (d-e) uj* for frictionless 
packings scales both with p and z — 4. Dashed lines indicate 
power laws with exponents as indicated. For details, see text. 



small Az are intimately connected to these floppy modes 
0, 1^ F'o^' frictionless systems they also cause the 

shear modulus G to become much smaller than the bulk 
modulus X — in fact, GjK ~ Az [l[ilEl- 

Our second main result is that we have found numer- 
ically that for frictional systems the ratio GjK only de- 
pends on Az (and not on e.g., /i), and for small Az it 
also scales as Az. To calculate the moduli, we start from 
the dynamical matrix which relates forces and displace- 
ments. Calculating, in linear order, the global stress re- 
sulting from an imposed deformation, the bulk and shear 
moduli are deduced [ijj. Since for Hertzian forces the 
effective spring constants scale as p^^'^, we have divided 
out this trivial pressure dependence. The results of our 
calculations are shown in Fig. 4. As could be expected, 
the (rescaled) bulk modulus K remains essentially con- 
stant. Surprisingly, the shear modulus G becomes much 
smaller than K for small p and large /i, and when plot- 
ted as function of Az, the ratio G/K \s found to scale 
as Az, as was predicted in lll|. Hence, in packings of 
deformable spheres, both uo* and G/K scale with Az, 
regardless of the presence of friction. 

Discussion — The sudden change in the DOS when 
increasing p from zero hints at the singular nature of the 
/i — > limit. On the one hand, the nature of the dy- 
namical matrix suddenly changes in this limit because 
the rotational degrees of freedom which are irrelevant for 
/i = turn on as soon as ^ ^ 0. On the other hand, 
we have recently found that the slower the packings are 
allowed to equilibrate during their preparation, the more 
the density of fully mobilized contacts, i.e., those for 
which — fiF'^, tends to increase; in fact especially for 
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FIG. 4; Scaling of bulk modulus K and shear modulus G as 
function of p and /i — similarly to the data for ui* , the trivial 
p^^^ dependence has been divided out. (a) The rescaled bulk 
modulus K (curve) essentially levels off for small p, while the 
shear modulus G (symbols as in Fig. 3) varies strongly with 
both p and ^. (b) G/ K scales like the excess contact number 
for small Az. 



isostatic limit. But there is an important difference to the 
frictionless case: while there the isostatic limit is auto- 
matically reached in the hard particle/small p limit, this 
is not necessarily so for the frictional case — here p and 
z are not directly related, and only for large friction does 
z approach isostaticity at small pressures. This isostatic 
point is relevant in practice: most materials have a value 
of of order 1, and as Fig. Eland 4 illustrate, one then 
observes approximate scaling over quite some range. 

We are grateful to M. Depken, L. Silbert, S. Nagel, 
D. Frenkel, H. van der Vorst and T. Witten for illumi- 
nating discussions. ES acknowledges support from the 
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foundation FOM, and MvH support from NWO/VIDI. 



small the fraction of fully mobilized contacts in slowly 
equilibrated samples becomes very substantial I l9l l20l |. 
The effect of these fully mobilized contacts on the DOS 
depends on additional physical assumptions. For exam- 
ple, if we assume that, for some reason, the contacts re- 
main constrained at the Coulomb threshold in the vibra- 
tional dynamics, we expect them to have an enhanced 
DOS for small pressures at all ji. More likely, these con- 
tacts would slip, leading to an initial strongly nonlinear 
response after which no contacts would be fully mobilized 
anymore, and our results for the DOS would go through 
essentially unchanged. 

Note that even a small difference between dynamic 
and static friction could suppress the effect of the fully 
mobilized contacts. Moreover, for realistic values of the 
friction (/i ^ 0.7, say) these effects are not very impor- 
tant since there the fraction of fully mobilized contacts 
is small. Thus, the results of this paper will apply di- 
rectly to packings with experimentally relevant values of 
the friction. 

A second issue that deserves further attention is the 
nature of the soft modes. Our scahng result for G/K sug- 
gest that for frictional systems these are dominated by 
shear-like (volume conserving) deformations, just as for 
frictionless systems. Apparently, rotations and particle 
motions couple such as to allow large scale floppy-mode- 
like distortions of frictional isostatic packings. Indeed, 
numerically obtained low frequency eigenmodes of fric- 
tional and frictionless systems look remarkably similar. 
Whether the scaling of w*, the scaling of G/K and the 
nature of floppy modes are similarly related in more gen- 
eral systems, such as packings of frictionless or frictional 
non-spherical particles, is an important question. 

Outlook — Our study of the density of vibrational 
states for frictional systems gives strong evidence for a 
scenario partly analogous to the one for frictionless pack- 
ings: frictional granular media become critical and ex- 
hibit scaling when their contact number approaches the 
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